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Résumé 

We propose a A. G. M. algorithm for the détermination of the characteristic poly- 
nomial of an ordinary non hyperelliptic curve of genus 3 over F 2 w . 

Nous nous proposons dans la présente note de donner un algorithme de calcul du 
polynôme caractéristique de la jacobienne d'une courbe de genre 3 ordinaire et non- 
hyperelliptique sur k = ¥ 2 n . Cet algorithme appartient à la catégorie des méthodes dites 
A. G. M. pour le calcul du nombre de points introduites par J.-F. Mestre dans les cas 
hyperelliptiques (cf. [2] et [3]). Nous illustrons son fonctionnement avec le calcul réalisé 
sous MAGMA en 25 heures pour une courbe sur F 2 72 . Une implémentation plus efficace 
(utilisant les méthodes de multiplication rapide) est en cours. Les démonstrations des 
résultats utilisés seront disponibles dans la thèse de l'auteur. Notons enfin que la méthode 
analogue dans le cas réel permet sans doute le calcul d'une matrice de Riemann de la 
jacobienne d'une courbe réelle, imitant en cela l'A. G. M. classique de Gauss dans le cas 
des courbes de genre 1. 

1 Méthode A. G. M. 

Etant donnée une courbe C de genre g sur k = ¥ 2 n ordinaire, la mise en œuvre d'une 
méthode A. G. M. pour la détermination du polynôme caractéristique procède globalement 
de la manière suivante : 

1. on relève la courbe sur une extension K non ramifiée de (Q>2- Ce relèvement ne doit 
pas être quelconque si l'on veut pouvoir effectuer les calculs dans une extension non 
ramifiée. Nous donnerons dans le cas du genre 3 au paragraphe 3, un bon modèle 
sur une extension de degré au plus AN de Q2. 

2. Une fois en caractéristique 0, par analogie avec le cas complexe, on calcule alors 
en fonction des coefficients de ce modèle 2 9 rapports de «thêta constantes» de la 
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forme suivante 
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dans le cas du genre 3 les 8 rapports qui nous intérèssent parmi les 2 9 ~ 1 (2 9 + 
1) = 36 non nulles sont exactement ceux pour lesquelles le rapport précédent est 
congru à 1 modulo 8. Ces derniers orientent en effet le noyau de l'isogénie pour que 
celle-ci coincide avec celui du relèvement du Verschiebung qui est un sous-espace 
totalement isotrope pour le couplage de Weil. Quitte à faire alors un changement 
de base symplectique on peut supposer que ceux-ci sont tels que / = 0. On note 
ces rapports (A^) e( z(z/2Z)s ■ 

Il faut exprimer ces rapports en fonction des coefficients de la courbe. Dans le cas 
hyperelliptique, les relations sont données par la formule de Thomae [4, livre II]. 
L'objet du paragraphe 2 est de donner un analogue de ces formules pour le cas 
du genre 3 non hyperelliptique. Alors que dans le cas hyperelliptique sont utilisés 
des invariants binaires reliés aux points de Weierstrass, on utilise ici des invariants 
ternaires basés sur des produits de déterminants de bitangentes. 
3. On duplique ces constantes par une généralisation de la moyenne arithmético- 
géométrique pour le genre 1 : 



la racine carrée de x G 1 + 80 étant choisie comme celle congrue à 1 modulo 4. 



-Ki sont les g racines du Frobenius inversibles modulo 2 ([3]). Ce théorème peut être 
démontré par la théorie de la multiplication complexe, l'argument de convergence 
provenant d'un théorème de R. Caris (cf. [1]). 

5. On calcule alors le polynôme minimal de (3 = a + 2 gN /a. Dans le cas g = 1,2 
on peut montrer que le polynôme minimal est essentiellement donné par ce nom- 
bre unique qu'il suffit de connaître avec une précision N/2 respectivement N ce 
qui implique respectivement 3iV*/2 et 2N itérations. Ce n'est plus le cas en genre 
3 : on utilise alors un analogue 2-adique de la méthode LLL pour la détermina- 
tion du polynôme minimal. L'utilisation brutale de cette méthode semble toutefois 
nécessiter la détermination de a avec une précision égale à UN. 

6. Lorsque g < 3 on montre alors que la connaissance de ce polynôme minimal permet 
de déterminer le polynôme caractéristique au signe près. Ce n'est plus toujours vrai 
en genre supérieur (cf. [3]). 

7. Il reste alors à lever l'ambiguité sur le signe. Cela peut se faire en effectuant un 
calcul dans la jacobienne (ce qui permet également de vérifier l'algorithme) ou si 
l'on connaît des points rationnels d'ordre 4. 





converge vers a = ±7Ti . . . ir g où les 
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Les deux paragraphes suivants donnent dans le cas g = 3 non hyperelliptique les 
principaux éléments qui permettent de mettre en œuvre la méthode décrite, en particulier 
la détermination des thêta constantes et le choix d'un bon modèle. Comme il a été dit, 
les justifications théoriques seront données ultérieurement. 



2 Détermination des thêta constantes 



Soit K est un corps algébriquement clos de caractéristique et C/K une courbe 
de genre 3 non hyperelliptique que l'on considère plongée dans P 2 comme une quartique 
lisse. Nous utilisons ici des travaux de Riemann et Weber (cf. [5] et [6]). Ces derniers sont 
basés sur le calcul des 28 bitangentes à C (i.e les droites l telles que (l • C) = 2P + 2Q) 
que nous explicitons comme suit. 

On note S = {L G Pic 2 (C),L 2 = JC} l'ensemble des fibrés thêta caractéristiques (où 
/C désigne le fibré canonique) qui se partagent en So et Si sous-ensembles pour lesquels 
h (L) = (resp. 1) et dits fibrés pairs (resp. impairs). On montre alors que l'ensemble 
des bitangentes est canoniquement en bijection avec Si. 

D'autre part, le choix d'une base symplectique de Jac(C)[2] permet d'associer à chaque 



point de 2-torsion une matrice 2x3 notée 



à coefficients dans F 2 (appelée car- 



actéristique) et de réaliser une bijection explicite avec S. On montre alors qu'à chaque 
bitangente est associée une caractéristique impaire (i.e. pour laquelle e ■ e' = 1). 
Afin de mieux comprendre la combinatoire de ces caractéristiques, on introduit un sous- 
ensemble de l'ensemble des caractéristiques appelé ensemble principal noté ([i])i=i ...7 
caractérisé par les propriétés suivantes : 

- Toute caractéristique impaire s'écrit soit [i] soit [i] + [7], i 7^ j. 

- Toute caractéristique paire s'écrit soit [0] soit [i] + [j] + [k] avec i,j, k distincts. 



On a par exemple 
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[4] = 



1 1 
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On note alors la bitangente associée à [i] et celle associée à [i] + [j] . 
Quitte à réaliser une transformation linéaire on peut supposer 



< 



/3i : xi = 
(3 2 : x 2 = 
: x 3 = 
P 4 : xi + x 2 + 23 = 



P 5 = a\x\ + a 2 x 2 + 0,3X3 
Pe = a'i^i + a' 2 X2 + 0-3X3 
/3r = a"xi + a 2 'x 2 + 0.3X3 
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on a alors le résultat suivant dû à Riemann [5] 

Théorème 1 II existe un plongement de C dans P 2 de la forme 

\/xmi + \/x 2 U 2 + V x 3^3 = 

où ui,U2,v>3 sont linéaires en xi,x 2 ,x 3 donnés par 

u± + u 2 + u 3 + xi + x 2 + x 3 = 
ÏÏT + « + ïïf + fca i x i + fc °2*2 + fca 3 x 3 = 

t + ? + ? + k '^ Xl + fc ' a 2 X2 + fc ' a 3 X 3 = 
$ + 1 + f + fc "°ï œ l + fe " a 2 X 2 + fc"03»3 = 

awec k,k',k" solutions des deux systèmes linéaires suivants : 
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5o«s cette forme une expression des 28 bitangentes est alors donnée par 



01 ■ Xi = : x 2 = /? 3 : x 3 = 
023 ■ ui = /?i 3 : u 2 = 0i2 : u 3 = 



/?4 : xi + x 2 + x 3 = /3 5 : a\X\ + a 2 x 2 + a 3 x 3 = 
06 : a[xi + a' 2 X2 + a' 3 x 3 = 0j : a'/xi + a 2 x 2 + a 3 x 3 = 

/3u : ni + x 2 + x 3 = 0i 5 : ^ + fca 2 x 2 + fca 3 x 3 = 



ai 

«1 







/3i 6 : ^ + fc'a 2 x 2 + fc'a 3 x 3 = /?i 7 : ^ + k"a' 2 f x 2 + ^''ag'xg = 

02a : xi + u 2 + x 3 = 
/3 26 : fc'aixi + ^ + À;'a 3 x 3 = 



025 : fcaiXi + ff + ^«3^3 = 

<J'„'I„. i «2 i !.//„//„ 



/3 3 4 : xi + x 2 + u 3 = 
/3 36 : k'a' lXl + A;V 2 x 2 + ^ = 



027 : fc'V/xi + ff + fc"a 3 'x 3 = 





/3 35 : fcaixi + ka 2 X2 + 
/3 37 : k"a'{ Xl + Fa'^ + ^ = 
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a^l-fc^) + a'J.d-W'a'ia'D + ^'(l-fc»^') = 
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Si on se donne C sous une forme y/xiu\ + y/x 2 u 2 + y/x 3 u 3 = il nous reste à montrer 
comment déterminer les On peut réaliser cela au moyen du calcul suivant 

1. On calcule Di(X), déterminant de la hessienne de la famille de coniques suivante 

Qi(X) = X 2 (x 2 u 3 ) + X(x 1 u 1 - x 2 u 2 - x 3 u 3 ) + (x 3 u 2 ). 

2. On calcule Ri(xi, x 2 , x 3 ) = Res(Di, Qi, À). 

3. On calcule de même R 2 {x\, x 2 , x 3 ) avec la famille 

Q 2 (X) = X 2 (x 1 u 3 ) + X(x 2 u 2 - xiui - x 3 u 3 ) + (x 3 ui). 

4. Enfin R = pgcd(i?i, R 2 )/(x 3 u 3 ) = ft 7 i=4 fc. 

La détermination des bitangentes est essentielle au calcul des thêta constantes. On a 
en effet le résultat suivant essentiellement contenu dans Weber (cf. [6]) 

Théorème 2 Soit [x] = [i] + [j] + [k] . Alors 

où l'on note [/3 h , (3 h , (3 h ] = det(/3 h , (5 h , (3 h ). 

3 Bon modèle 

On reprend ici les notations du paragraphe 1. Une courbe de genre 3 non-hyperelliptique 
sur k = ¥ 2 n a au plus 7 bitangentes et le modèle de Riemann ^/x\û[ + y/x 2 u 2 + y/x 3 u 3 = 
n'est plus valable en caractéristique 2. On commence donc par chercher une description 
simple des courbes qui nous intéressent 

Lemme 1 Toute courbe de genre 3 ordinaire, non hyperelliptique sur k est isomorphe à 
une courbe d'équation 

(ax 2 + by 2 + cz 2 + dxy + exz + fyz) 2 — xyz(x + y + z) = (1) 

avec la condition suivante 

abc(a + b + d)(a + c + e)(b + c + f)(a + b + c + d + e + / + 1) / 0. 

Inversement toutes les courbes qui vérifient ces conditions sont des courbes de genre 3 
ordinaires et non hyperelliptiques. 

De plus l'isomorphisme est défini sur une extension au plus cubique de k. 

On considère donc C donnée par les coefficients a, b, c, d, e, f ci-dessus. Nous voulons 
remonter ce modèle sur une extension de Q 2 non ramifiée de telle sorte que ses bitan- 
gentes soient faciles à déterminer, c'est-à-dire qu'elles soient définies par les formules du 
théorème 1. Nous allons pour cela légèrement déformer le modèle de Riemann. 
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Tout comme la forme y 2 = x 3 + ax + b doit être transformée en caractéristique 2 en 



y +xy 



■ax + b, nous allons transformer notre modèle de Riemann en rajoutant des 



termes à la Artin-Schreier. Remarquons pour cela que le modèle de Riemann ^x\û[ + 
yJX2U2 + \Jx 3 u 3 = s'écrit également 

>i 2 



X\U\ 

X2U2 
X3U3 

Y 1 + Y 2 + Y 3 
Considérerons le modèle suivant sur k 



Y 2 

ï 2 

Y 2 



Y? 

'2 



hYi = hvi 

Y'4 + hY2 = hv 2 
Y 2 + l 3 Y 3 = l 3 v 3 

h + l2 + h = 
[Y 1 + Y 2 + Y 3 = l 

où h,h,l2,vi,V2,v 3 et l sont linéaires. On a par exemple 

Proposition 1 C donnée par le modèle (1) est isomorphe sur k à la courbe donnée par 
le modèle ci-dessus avec 

h = x,l 2 = y,k = x + y 
l = z 

vi = bcy + (c + f)z 
V2 = acx + dey + (c + f)z 
{v 3 = acx + (d + b)cy + (1 + e + c + f)z 

On relève alors les expressions obtenues sur K à l'identique sauf pour £3 = — 2z — x — y 
et on écrit 

\2Y l + l 1 ) 2 = l l (Av l + l 1 ) 
(2Y 2 + l 2 ) 2 = l 2 {4v 2 + l 2 ) 
{2Y 3 + l 3 ) 2 = l 3 {Av 3 + l 3 ) 
(2Yi + m) + {2Y 2 + u 2 ) + {2Y 3 + u 3 ) = 

xi,y = x 2 ,z = -(xi + X2 + x 3 )/2 on a 



Si on effectue le changement de coordonnées x 
le modèle de Riemann 



C : ^xi^ivi + h) + y/x 2 (4v 2 + h) + ^Jx 3 {Av 3 + l 3 ) = 

4 Exemple 



(2) 
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On considère la courbe suivante : 



C : (x 2 + uj 2 y 2 + oj^z 2 + ojxy + (uj 2 + l)xz + (uj d + l)yz) 2 — xyz(x + y + z) 
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définie sur k = F 2 n, N = 72 où lo engendre le groupe multiplicatif k*. 

Les calculs sont réalisés grâce au logiciel MAGMA version 2.9 sur un Pentium III à 1.13 

Ghz avec 2 GigaOctets de mémoire centrale. 

On calcule alors les vi,V2,vs du modèle de la proposition 1 

Vl = uj 6 y + (lo 4 + lo 3 + l)z 

v 2 = lo 4 x + Lo 5 y + (lo 4 + lu 2 + l)z 

v 3 = oj 4 x + (lo 6 + uj 5 )y + (uj 4 + lo 3 + uj 2 + l)z 

On note K l'extension de degré d non ramifiée de Q2, O son anneau d'entiers et vu £ O 
qui se réduit sur lo. Le modèle (2) de C est défini par 

Ul = (-2w 4 - 2w 3 - l)xi + (4w 6 - 2w 4 - 2w 3 - 2)x 2 + (-2w 4 - 2w 3 - 2)x 3 
u 2 = (2w 4 - 2w 2 - 2)xi + (4w 5 - 2w 4 - 2w 2 - l)x 2 + (-2w 4 - 2w 2 - 2)x 3 
u 3 = (2w 4 - 2w 3 - 2w 2 - 2)xi + (4u> 6 + 4w 5 - 2w 4 - 2w 3 - 2w 2 - 2)x 2 + 
(-2w 4 - 2w 3 - 2w 2 - l)x 3 

On effectue le calcul des 4 bitangentes qui nous manquent grâce à l'algorithme du para- 
graphe 2 (temps : 1/2 heure, elles sont définies sur K dans notre cas) puis de l'ensemble 
des bitangentes et de tous les rapports non nuls 
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grâce aux théorèmes 1 et 2. On sélectionne les 8 rapports congrus à 1 modulo 16 et on 
calcule leur racine carrée dans K (temps : 15 mn) 

'l + 0(2 4 ) 

1 + (w 16 + w 14 + w 12 + w 11 + w 9 + w 5 + w 4 ) ■ 2 3 + 0(2 4 ) 
1 + {w 16 + w 14 + w 12 + w 11 + w 10 + w 8 + w 4 + w 3 ) ■ 2 3 + 0(2 4 ) 
1 + (w 10 + w 9 + w 8 + w 5 + w 3 ) ■ 2 3 + 0(2 4 ) 
' 1 + (w 6 + w 4 + w 3 ) • 2 3 + 0(2 4 ) 
1 + (w w + w 14 + w 12 + w 11 + w 9 + w 6 + w 5 + w 3 ) ■ 2 3 + 0(2 4 ) 
1 + (w 16 + w 14 + w 12 + w 11 + w 10 + w 8 + w 6 ) • 2 3 + 0(2 4 ) 
1 + (w 10 + w 9 + w 8 + w 6 + w 5 + w 4 ) • 2 3 + 0(2 4 ) 

On effectue alors l'itération de la formule de duplication 12./V fois (temps : 24 heures) 
qui nous donne au signe près la valeur du produit des racines inversibles du Frobenius 
modulo 2 

= i + 28 + 2 9 + 2 11 +2 13 + 2 15 + 2 16 + 2 17 + 2 18 + 2 21 +2 23 + 2 24 +. . . + 2 787 + 2 790 + 2 791 +O(2 793 ) 



7 



On trouve alors grâce à LLL (temps : 1 seconde) le polynôme minimal de (3 



P SJ m{x) = x 4 - 52767044410803560460262696266497a; 3 - 

78121277277710794719527572033891108646286909 • 2 7 V + 
610161746623391968394415142270976679056928051874538813 • 2 2 ' 72 x + 
46918330565326150855288775851644884890720289023905899509851903489 • 2 3 ' 

D'où le polynôme caractéristique au signe près 

P{x) = x 6 - 9925657555a; 5 + 1108548370771462406931x 4 

-146512229527151304651245280013057x 3 + 1108548370771462406931 • 2 7 V 
-9925657555 • 2 2 ' 72 x + 2 3 ' 72 . 

Pour déterminer le signe, on étudie alors les points d'ordre 4 de la jacobienne. Plus 
précisément dans notre cas, 2 6 divise -P(l) mais 16 ne divise pas P(— 1). Il suffit donc 
que la jacobienne de C possède un point d'ordre 4 rationnel pour conclure que P(X) est 
le «bon» polynôme. Considérons par exemple le point d'ordre 2 (P± + Q\ — (P2 + Q2)) = 
^((x-C)-(y-C)). On cherche un point de la jacobienne sous la forme D = P+Q+R—3P2 
tel que 2D = P 1 + Q 1 - (P 2 + Q 2 ) ou encore 2{P + Q + R)- (P 1 + Q 1 + 5P 2 - Q2) = 0. En 
utilisant par exemple MAGMA, on montre que l'espace L(P\ + Q\ + 5P 2 — Q2) qui est 
de dimension 4, contient une fonction définie sur F 2 72 ayant trois zéros d'ordre 2, d'où 
l'existence d'un point d'ordre 4 rationnel sur F 2 72 (temps : 10 secondes). 

Remarque : A posteriori on constate que la précision nécessaire est seulement de 
10-72, le calcul est alors réalisé en 21 heures. Il serait donc très utile d'avoir une meilleure 
borne pour LLL. 

Références 

[1] R. Caris : Approximation of canonical lifts, in préparation, 2002 disponible sur 
http ://www. math. leidenuniv.nl/ caris/. 

[2] J.-F. Mestre : lettre à Gaudry et Harley, disponible sur 
http ://www. institut. math.jussieu.fr/, (2000). 

[3] J.-F. Mestre : Algorithmes pour compter des points en petite caractéristique en 
genre 1 et 2, disponible sur www.maths.univ-rennesl.fr/crypto/2001-02/mestre.ps. 

[4] D. Mumford : Tata Lectures on Thêta, Vol 1,2, Birkhâuser, (1983). 

[5] B. Riemann : sur la théorie des fonctions abéliennes, Oeuvres de Riemann, Paris : 
Gauthier- Villars, (1898), deuxième édition, p. 487. 

[6] H. Weber : Théorie der abelschen Functionen vom Geschlecht 3, (1876). 



8 



